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SUMMARY 

An analysis is made for laminar forced-convection heat transfer on 
a flat plate with unsteady surface temperature. The deviation of the 
instantaneous heat-transfer rate from the quasi-steady value is computed. 
Results are obtained for Prandtl numbers in the range 0.01 to 10. The 
deviations from quasi- steady heat transfer increase markedly with in- 
creasing Prandtl number. The findings reported here should apply approx- 
imately in the entrance region of ducts and should also provide an upper 
bound on deviations from turbulent quasi-steady heat transfer. 


INTRODUCTION 

It is often necessary to compute the forced-convection heat trans- 
fer from a surface whose temperature is changing with time. This problem 
is much simplified when it is supposed that the boundary layer passes 
through a succession of instantaneous steady states. Such a boundary 
layer is called quasi- steady . Under this assumption, the heat transfer 
is computed by instantaneous application of steady-state heat-transfer 
relations. The quasi-steady supposition is also invoked in heat-transfer 
experiments employing the transient technique, where the instantaneous 
measurements are used to determine steady-state coefficients. 

In reality, there is always a difference between the actual instan- 
taneous heat transfer and the quasi-steady value. The extent of the 
deviation depends upon the response characteristics of the boundary 
layer, as well as on the rapidity of the changes in surface temperature. 

The aim of this analysis is to find the first- and second-order de- 
viations of the actual instantaneous heat transfer from the quasi-steady 
value and to learn how these deviations depend upon the Prandtl number. 
The final results provide a rapid and accurate quantitative means for 
determining when a given set of surface temperature and free-stream ve- 
locity data lead to essentially quasi-steady heat transfer. 
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The system chosen for study is a flat plate alined parallel to a 
steady laminar flow as pictured in the following sketch: 


v 



The surface temperature T w is spatially uniform but is permitted to 
take on arbitrary, but continuously differientable, variations with time. 
The free-stream temperature T„ is taken to be constant. Results are 
given here for fluids having Prandtl numbers between 0.01 and 10. This 
investigation constitutes an extension of previous work for Pr = 0.72 
(air and other gases) reported in reference 1. Those primarily inter- 
ested in results are invited to pass over the section on ANALYSIS. 

Readers interested in non- quasi -steady boundary layers are referred 
to the work of Moore and Ostrach (refs. 2 to 5), who studied the effects 
associated with time variations in free-stream velocity. 
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SYMBOLS 

Cp specific heat at constant pressure 

F Blasius velocity function 

h local heat- transfer coefficient, q/AT 

k thermal conductivity 

M free-stream Mach number 

ao 

Pr Prandtl number, c p |J./k = v/a 

q local heat-transfer rate per unit area 

R recovery factor 

T static temperature 

AT temperature difference j T w - T aw with frictional heating, 

T w - T„ without frictional heating 

• • • 

T W ,T W time derivatives of wall temperature 
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t 

Uoo 

u 

V 

X,x 

Y 


time 

free- stream velocity 

velocity component in x-direction 

velocity component in y-direction 

coordinate measuring distance along plate from leading edge 


transformed coordinate. 



y 

a 


coordinate measuring distance normal to plate 
thermal diffusivity, k/pCp 


expansion parameters, 


M T v - T - 


Y* ratio of specific heats 

Y 

T[ Blasius similarity variable, — 

6 dimensionless temperature, (T - T m )/(T y . - 1 ^) 

©1,02 functions of rj 
\ i absolute viscosity 

v kinematic viscosity 

p density 

t time 

$ function of t) 

Y stream function 

Subscripts : 

aw adiabatic wall 



inst 


instantaneous 
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qs 

quasi- steady 

w 

wall 

00 

free stream 


ANALYSIS 

Governing Equations 

The analysis begins with the equations expressing conservation of 
mass, momentum, and energy for unsteady laminar boundary- layer flow over 
a flat plate: 


a 


(p u ) + ^ (pv) + 


Bp _ 
St - 


0 



The boundary conditions appropriate to the problem are 

'S 

u -* U« 

> y = o 


u = o 
v = o 
T = T w (t) 


T 


y 


(1) 

( 2 ) 

( 3 ) 

( 4 ) 


To simplify the treatment, the variation of fluid properties in 
liquids will be neglected, while for gases it will be supposed that the 
properties may be approximated by pp = constant, pk = constant, Cp = 
constant. Fluid-property variations are not expected to have a decisive 
effect on the final results (which are expressed in the form of ratios). 


The equation for conservation of mass is satisfied by a stream func- 
tion Y defined as follows (ref. 2, eq. (5)): 


. R. By 
p 3y’ 
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Then, by replacing u and v in favor of Y and introducing the fol- 
lowing new variables 


« 



X = x, 



equations (2) and (3) may be rephrased as 


t = t (5) 


y yt + Wy ~ ¥ x ¥ yy = v °° ¥ yyy 


+ ^ ip 


w 


L w V ( 

w i y n V G — -J 

- T_ Y X " X Y " PrTYY r T - T 




, \2 

~ ^ YY' 


w 


(2a) 

(3a) 


For brevity, the derivatives with respect to X, Y, and t have been 
denoted by subscripts. In terms of the new variables, the boundary con- 
ditions (4) become 



Y 

9 



(4a) 


Since the transformation has removed all thermal effects from the momen- 
tum equation (2a), it is clear that y Y t - 


Solutions 

As has already been noted, the object of this analysis is to in- 
vestigate the deviations of the actual instantaneous heat transfer from 
the quasi-steady value. With this in mind, it is natural to seek a solu- 
tion for the temperature distribution in the form of a series expansion 
about the quasi-steady state. In reference 1, the following quantities 
were introduced to serve as expansion parameters for the series solution: 



Physically, these parameters represent the ratio of the time required 
for changes in surface temperature to diffuse across the boundary layer 1 

-'-The time required for changes to diffuse across a layer of thick- 
ness A is proportional to /oc, 1 wh i le the thermal boundary— layer thick- 
ness A is proportional to V m / U « • Hence, the diffusion time is pro- 
portional to x/U^. 
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to a time that is characteristic of the rapidity of the surface temper- 
ature changes. Thus, the p n serve as a measure of the promptness with 
which the boundary layer responds to impressed variations of surface 
temperature . 

Using the P n , the series solution for the dimensionless tempera- 
ture 0 is written as 

U£ 

9 = ©gCil) + Pl^l^) + 02^2^) + • • • + 2c (t - T — 5 < ^ T *^ ( 7 ) 

where is the well-known Blasius similarity variable given by 



The functions 9q and $ are associated with the quasi-steady tempera- 
ture distribution. When the Pn are small (corresponding to prompt 
response to impressed changes), the state is essentially quasi-steady. 

The stream function I is written in -terms of the Blasius variable 

F as 

(9) 

When the expressions for 6 and Y are substituted back into equa- 
tions (2a) and (3a) , and after terms are grouped according to whether 
they are multiplied by P^> Pg* • • • ’ following se l ordinary dif- 

ferential equations results: 



F "’ + 

FF" = 

0 


F(0) = F'(0) . 

= 0, F'(«) = 2 

( 10 ) 


9 o + PrF0^ 

= 0 

0 O (°) = 

H 

CD 

O 

IP 

II 

O 

(11) 


+ Prj 

CD 

H- 

i 

2F' 

CD~ 

i — 1 

CD 

)-° 

0l(O) = 0 1 (oo) = 0 

(12) 

®2 

+ Pr| 

- 

4F ' 

'0 2 - 40- 

) = o 

© 2 ( 0 ) = 02(°°) = 0 

(13) 


<t>" + 

PiTpP 

' + 

| (F") 2 " 1 

l-° 

4>(o) = $(-) = 0 

(14) 


A solution of equation (10) was first obtained by Blasius in 1908, 
was necessary to re- solve this equation to greater accuracy for 
present purposes. Numerical solutions to equations (ll), (l2), and (13) 
have been carried out as part of this investigation on an IBM 650 
Magnetic Drum Data-Processing Machine for Prandtl numbers of 0.01, 0.72, 
1.0, and 10. The function $ is associated with the aerodynamic heating. 


NACA TN 4311 


7 


For laminar conditions, this effect of viscous dissipation is usually of 
importance only for high-speed gas flows. Solutions of the aerodynamic- 
heating problem for the Prandtl number range of gases are available in 
the literature (e.g., ref. 6). 

These numerical results will be utilized in the heat- transfer cal- 
culations that follow. 


HEAT-TRANSFER RESULTS 


The instantaneous local heat flux at the plate surface q ins t 
be calculated by applying Fourier ! s law: 


■ ■ (> 4 


y=o 

After introducing the series expansion of (7) and taking account of the 
transformed variables of equations (5) and (8), the expression for q 
becomes 


linst - 


! $ < T » - T - 


®q(o) + ^(o) + V 2 e' z (o) + . . . + 


u. 

2 = p ( t , - t j 


4** (0) 


(15) 


where 0 q(o), 0j(O), . . . are abbreviations for (d^o/dt]) q. 


For liquids, the effects of aerodynamic heating will be neglected, 
and consequently the last term on the right side of the equation is omit- 
ted. For gases, it is convenient to introduce the adiabatic wall tem- 
perature T aw by the relation 


T = T 
x aw - L o 


1 + R 


(r ~ ll m 2 = 


me] = Too(l 


U* 


+ R 


2c T 


( 16 ) 


where R, the recovery factor, has been given in reference 1 as 

>*(0) 

Koy 


R- - 14 °) 
o v 


(16a) 
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Making use of equations (l6a), (l6), and (6), it is easy to show that 
the final expression for !i ns t ^ as same for liquids and for 

gases, namely, 


^inst 





e^(o) + . 


where 


AT = T - T 

i w aw 


for gases 
for liquids 

The quasi-steady heat transfer q qs is given by 


AT = T v - T< 


d qs 


■ -mP 

2 V vx 


AT G' 0 { 0) 


(17) 


(18) 


The important relation between the instantaneous and the quasi- 
steady heat transfer is then found by combining equations (17) and (18) : 



The quantities 6^(0), 0|(O), and e^ 0 ) have been found as P art of the 
solution of equations (ll), (12), and (13) and are listed in table I. 

For convenience, the ratios 0 ^(o)/0q(o) and 9 '^{o) / 9 '^(0) are plotted in 

figure 1 as a function of Prandtl number. It is seen from equation (19) 
that the deviations from quasi- steady heat transfer depend directly on 
the magnitudes of these ratios. Figure 1 shows that both ratios in- 
crease markedly with Prandtl number, demonstrating, for example, that 
liquid metals (Pr = 0.01) are less likely to experience deviations from 
quasi-steady heat transfer than is water (Pr * 5) . This trend might have 
been intuitively expected on physical grounds as a consequence of the 
relatively high thermal diffusivity of liquid metals. 


As alternative form of the results may be obtained by introducing 
heat-transfer coefficients as follows: 




qs 


= 

AT 


v _ q inst 
n inst ’ 


( 20 ) 
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Then, equation (19) may be rephrased as 



Equation (2l) is in a form useful for interpretation of heat-transfer 
coefficients obtained under transient conditions. 

To further facilitate calculations, a plot of the quasi-steady heat- 
transfer relation (eq. (18)) is presented in figure 2. 


Criterion for Quasi-Steady Heat Transfer 

From the series nature of the solution, it is expected that equa- 
tion (19) would be most accurate for small deviations of q.i ns t/cLqs from 
unity. This suggests that this equation can serve as an accurate and 
rapid means for checking whether a given situation may be treated as 
quasi- steady. When the free-stream velocity and surface temperature data 
of a particular situation lead to linst/^-qs * then that situation can 
be taken as quasi-steady for heat-transfer purposes. 

In a great many forced- convection applications, the quotient x/U^ 
(in seconds) will be a small number. As a consequence, the second term 
in the brackets of equation (19) will often be very small compared with 
the first term. Under these circumstances, it is possible to establish 
a simple criterion for determining when heat transfer from a surface with 
time-dependent temperature can be computed with sufficient accuracy from 
quasi-steady relations. Suppose that it is decided that an accuracy of 
5 percent is sufficient for many applications. Then, the upper curve of 
figure 3 distinguishes the conditions under which deviations from the 
quasi-steady state may be ignored in the computation of local heat trans- 
fer. Alternatively, if an accuracy of 2 percent is required, the lower 
curve of figure 3 gives the conditions under which quasi-steady relations 
are adequate. As a simple example, suppose that water (Pr = 5) is flow- 
ing at 50 feet per second over a flat plate. If T ir - T as 100° F and 
x = 1 foot, then the departures from the quasi-steady state will have 
negligible effect (^:2 percent) on the local heat transfer when the rate 
of change of surface temperature does not exceed 23° F per second at 
x = 1 foot. 


CONCLUDING REMARKS 

Although the analysis given here is for laminar flow over a flat 
plate, the findings may have a wider utility. Since the response of a 
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turbulent flow should be more rapid than that of laminar flow, the pres- 
ent results should give an upper bound on the deviations from turbulent 
quasi-steady heat transfer. Further, it is expected that the results 
reported here will apply approximately in the entrance region of ducts. 

It is felt that the final results should not be strongly sensitive 
to fluid-property variations, although caution should be exercised when 
fluids with unusual property variations are involved (e.g., near the 
critical point) . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 28, 1958 
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TABLE I. - FUNCTIONS NEEDED IN HEAT- 
TRANSFER COMPUTATIONS 


Pr 

0^(0) 

©2.(0) 

0 ' 2 ( 0 ) 

0{(O) 

9 2< 0) 

S^ol 

e^oT 

O.Ol 

-0.1032 

-0.1169 

0.02034 

1.133 

-0.1971 

.72 

-.5913 

-1.416 

.4739 

2.395 

-.8015 

1.0 

-.6641 

-1.751 

.6453 

2.637 

-.9717 

10 

-1.456 

-7.993 

6.140 

5.489 

-4.216 


(0)/9q(0) 
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Figure 2 


Quasi -steady heat -transfer results 
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NACA - Langley Field, V*. 



Figure 3. - Values of 
eq. (19)). 
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